#!/usr/bin/env python3
# -*- coding: utf-8 -*-
"""
Created on Sun Jun 30 21:35:37 2024

@author: liqingsimac
"""
import numpy as np
from scipy.integrate import odeint

def rhs(z,t):
    return [ z[1], (1-z[0]**2)*z[1]-z[0] ]
t=np.linspace(0,10,101)
y0=np.array([2,3])
y=odeint(rhs,y0,t)
import matplotlib.pyplot as plt
plt.plot(y[:,0],y[:,1],'.-')


